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Abstract. With a fantastic sensitivity improving significantly over the advanced GW detec¬ 
tors, Einstein Telescope (ET) will be able to observe hundreds of thousand inspiralling double 
compact objects per year. By virtue of gravitational lensing effect, intrinsically unobservable 
faint sources can be observed by ET due to the magnification by intervening galaxies. We 
explore the possibility of observing such faint sources amplified by strong gravitational lens¬ 
ing. Following our previous work, we use the merger rates of DCO (NS-NS,BH-NS,BH-BH 
systems) as calculated by Dominik et al.(2013). It turns out that tens to hundreds of such 
(lensed) extra events will be registered by ET. This will strongly broaden the ET’s distance 
reach for signals from such coalescences to the redshift range z = 2 — 8. However, with 
respect to the full inspiral event catalog this magnification bias is at the level of 0.001 and 
should not affect much cosmological inferences. 
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1 Introduction 

Gravitational waves (GW thereafter) are expected to be registered soon by the upgraded 
advanced LIGO/VIRGO detectors [1, 2]. This will open an era of experimental gravitational 
wave astrophysics with the expected detection rates of order ranging from tens to a thousand 
events per year. The real breakthrough will come with the new generation of detectors, an 
example of which — the Einstein Telescope (ET) has already went through the initial design 
study [3] . Because such an instrument will improve an order of magnitude in sensitivity over 
the advanced laser interferometer detectors LIGO and VIRGO, simple scaling arguments 
lead us to expect tens to hundreds thousands of detections per year. With such a number 
of detectable events and with the detector horizon reaching 1-2 Gpc one may expect 
that non-negligible number of sources would be gravitationally lensed by galaxies lying in 
between. We considered this problem in our previous papers [4, 5] . In particular the analysis 
performed in [5] was very comprehensive in the sense of taking into account full population 
of double compact objects (DCO) i.e. NS-NS, NS-BH, and BH-BH binaries, as well as taking 
the cosmological merger rates at different redshifts as suggested by the population synthesis 
model (using StarTrack code) in Dominik et al. [6]. The result was that ET would register 
about 50 - 100 strongly lensed inspiral events per year with statistics dominated by the BH- 
BH systems. These results suggest that ET will provide a considerable catalog of strongly 
lensed events during a few years of its successful operation. Our previous estimates in [5] 
were obtained under assumption that DGO systems intrinsically have signal to noise ratio 
(SNR) greater or equal to the threshold po = 8. This was a reasonable assumption since 
we wanted to estimate the rates of lensed events in the population of sources detectable to 
the ET. However, gravitational leasing effect will magnify the amplitudes of the GW sources 
increasing this way the SNR of each lensed DCO system. Therefore, intrinsically faint sources 
(having SNR < 8 ) may now become observable to the ET. In this paper, we supplement our 
previous study by considering sources which are intrinsically below the detection threshold 
and consequently would not be observed had not they been magnified by the lens. 

The paper is organized as follows. In Section 2, we briefly recapitulate our methodology 
(referring to [4, 5] for detailed calculations) and review DCO catalog build from evolutionary 
population synthesis code which we use thereafter. Then, in Section 4 we present our results 
and conclusions. 

Throughout the paper we comply with the notation and nomenclature of [5] and [6]. 
For the sake of consistency with previous works, we assume flat FRW cosmological model. 
In particular the expansion rate in this model reads: 


H{z) = + z)^ + (1 - nm) 


( 1 . 1 ) 
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with Hq = 70 km/s/Mpc, Q,m = 0.3 as in Dominik et al. [6]. We will also adopt the 
notation: E{z) = H{z)/Hq and f = fg — the non-dimensional comoving distance. 


2 Methodology 


Our goal is to estimate the number of GW sources which would be magnified above the 
detector’s threshold. However, unlike in our previous estimates we will admit them being 
intrinsically faint, i.e. such that they would not be detected haven’t they been lensed. From 
the physical point of view a relevant observable quantity is the dimensionless amplitude of 
the GW - the strain h{t) = F+/i+(f) -|- Fxhxit), where -|- and x denote two independent 
polarizations. However, the detector is not able to register pure signal, but rather the signal 
is buried in the detector’s noise, thus we have to use matched filtering technique in which 
signal to noise ratio (SNR) p is the integrated signal spectral power weighted down by noise 
power spectral density. In particular, for a single detector it reads (for more details see [4] 
and references therein): 



where: di is the luminosity distance to the source, 0 is the orientation factor capturing part of 
sensitivity pattern due to (usually non-optimal) random relative orientation of a DCO system 
with respect to the detector, tq above is detector’s characteristic distance parameter. After 
Taylor & Gair [7] we consider two options: the ET initial design, which gives tq = 1527 Mpc 
and the advanced “xylophone” configuration, which gives tq = 1918 Mpc. We also assume 
that Cifmax) = 1 (for justification and more details see [7]). The orientation factor 0 is 
defined as 

0 = 2[Fl{l + cos^ if + 4F2 cos^ if/^ (2.2) 

where F^ = ^(1 + cos^ d) cos 2(/)cos 2^/' — cos 0 sin 2(/)sin2'(/; and F^ = |(1 + 

cos^ 0) sin 2(;i)cos 2^ -|-cos 0 sin 2(^ cos 2^ are the interferometer strain responses to different 
polarizations of gravitational wave. 

Probability distribution for 0 calculated under assumption of uncorrelated orientation 
angles {6, (f>, -ijj, t) is known to be of the following form: 

PeiQ) = 50(4 - ef/256, if 0 < 0 < 4 (2.3) 

Pe{Q) = 0, otherwise 


The yearly detection rate of DGO sources originating at redshift Zs and producing the 
signal with SNR exceeding the detector’s threshold po (previously assumed that sources 
intrinsically have SNR parameter p = 8 but now we relax this assumption) can be expressed 
as: 


N{> pq\zs) = 


f 


dN{> po) 
dz 


dz 


(2.4) 


where 


dN{> po) 
dz. 


= 47r 


c 


nojzs) r'^izs) 

1 + Zs F{Zs) 


Ce{x{zs,po)) 


(2.5) 


is the rate at which we observe the inspiral DCO events (sources) that originate in the redshift 
interval [z, z + dz]. In Eq.(2.5) no(zs) denotes intrinsic coalescence rate in the local Universe 

at redshift 2^, Ceix) = Pe{&)dQ and x{z,p) = f (1 -|- ' 
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Here, like in the previous paper [5], we use the values of inspiral rates no(zs) reported 
by [6] for each redshift slice they considered. To be more specific they evolved binary systems 
from ZAMS until the compact binary formation under certain well motivated assumptions 
about star formation rate, galaxy mass distribution, stellar populations, their metallicities 
and galaxy metallicity evolution with redshift (“low end” and “high end” cases). In order to 
make straightforward comparison with the results presented in [5] we consider all evolutionary 
scenarios leading to the DCO formation, presented in Dominik et al. [6] that is: the standard 
scenario, optimistic common envelope scenario, delayed SN explosion and high BH kick 
scenario. We have taken the data from the website http:www.syntheticuniverse.org, 
more specifically the so called “rest frame rates” in cosmological scenario. For the chirp 
masses we have assumed the following values: 1.2 Mq for NS-NS, 3.2 Mq for BH-NS and 
6.7 Mq for BH-BH systems. According to Dominik et al. [8], they represent average chirp 
mass for each category of DCO simulated by population synthesis. 

Since we relax the fixed value of the intrinsic signal to noise ratio, instead of Eq. (2.5) 
we have to start with the differential inspiral rate per redshift and per SNR parameter p: 


d'^N 

dzgdp 



nojzs) r'^izs) 

I + Zs E{Zs) 


Pq{x{zs,p)) 


x{Zs,p) 

p 


( 2 . 6 ) 


Concerning gravitational leasing we adopt the same approach as in our previous paper 
[5], i.e. we assume conservatively that the population of lenses comprise only elliptical galax¬ 
ies. Therefore, we will model the lenses as singular isothermal spheres (SIS) which is a good 
approximation of early type galaxies [11]. Characteristic angular scale of lensing phenomenon 
is set by the Einstein radius, which for the SIS model reads: 9e = dvr (f)^ where a 

is the velocity dispersion of stars in lensing galaxy, dA{zi, Zg) and dA{zs) are angular diameter 
distances between the lens and the source and to the source, respectively. It is convenient 
to use the Einstein radius as a unit and convert the angular distance of the image (w.r.t. 
the center of the lens) 9 or the angular position of the source (5 to respective dimensionless 
parameters: x = y = Then the necessary condition for strong lensing (multiple 
images) is y < 1. Images (brighter /+ and fainter one /_ ) form at locations x± = 1 ± y 
with magnihcations: p.± = ^ ± 1. Hence, the gravitationally lensed GW signal would come 
from these two images with appropriate relative time delay (see [4, 5]) and with different 

amplitudes: h± = y/JI± h{t) = ± 1 h{t) where h{t) denotes the intrinsic amplitude (i.e. 

the one which would have been observed without lensing). Analogous relations are valid for 
the SNR parameter p. Assuming the threshold SNR for detection po = 8j one can observe 
lensed images (of the source with an intrinsic SNR equal to Pintr.) if the misalignment of the 
source with respect to the optical axis of the lens satishes: 

-1 

(2.7) 


y± — y±,max — 


Pintr. 


Tl 


These conditions (for the /+ and /_ image) influence elementary cross section for lensing 
(see e.g.[4]): 


Scr,±{^, zi,zs,p) = TT9%y'i „^^^ = levr^ 


7/2 

yzL^max 
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Figure 1. The observed lensed GW event number distribution as a function of p. “Low-end” 
metallicity galaxy evolution and standard model of DCO formation are assumed. 


which is necessary to calculate optical depth for leasing leading to magnifications of /+ and 
/_ images above the threshold: 


r±{zs,p) = 




da dvr ( 

ddo 


1^2 rjiii 

(2.9) 


In analogy to and in order to comply with our previous papers, [4] and [5], we model the 
velocity dispersion distribution in the population of leasing galaxies as a modified Schechter 

function ^ (^) (^) ) with the parameters n^,a^,a and /3 we taken 

after Choi, Park &: Vogeley [12] (for discussion about such choice in view of other data on 
velocity dispersion distribution functions see [5]). Optical depth for leasing depends on the 
survey duration. However, it turns out (e.g. our previous papers) that the survey time (i.e. 
1 year, 5 year and continuous) has little effect on the detection rate. Thus, in this work we 
only consider the case of continuous search. 

We can now combine Eq.(2.6) and Eq.(2.9) in order to calculate quantities like the 
leasing rate of intrinsically faint {p < po) DCOs having or /_ images magnified above the 
threshold pQ\ 

r ^max . rp° , , d'^N , 

Nlensed,±= / dZg / ^ dp (2.10) 

Jo Jo 


to JO dzsdp 

or differential leasing rates with respect to p or Zg respectively: 


dNi 


ensed,zb 


dp 

^^lensed,:L 


i 


r±(zg,p) 


d^N 


dzg 


dzg 


dzgdp 

r I c"'' j 


( 2 . 11 ) 


( 2 . 12 ) 


3 Results and discussion 

Table 1 shows the expected yearly detection rate of lensed DCO inspiral events having in¬ 
trinsic SNR below the threshold of detection pintr. < Po = 8 and magnified strongly enough 
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Table 1. Expected numbers of lensed GW events with pintr < 8 for which the /_ image is magnified 
above threshold po = 8. We also assumed the continuous survey, which in practice means that its 
duration is longer than 5 years. Nomenclature of DCO formation scenarios and galaxy metallicity 
evolution follows that of Dominik et al. [6] . Predictions for two configurations of the ET are given. 


DCO scenario 
metallicity evolution 

standard 
High; Low 

optimistic CE 
High; Low 

delayed SN 
High; Low 

high BH kicks 
High; Low 

NS-NS 
initial design 
xylophone 

0.6; 0.4 
0.9; 0.8 

5.1; 5.3 

9.4; 9.6 

0.6; 0.5 

1.0; 0.9 

0.6; 0.4 

1.0; 0.8 

BH-NS 
initial design 
xylophone 

1.2; 1.4 
1.4; 1.5 

2.2; 2.1 

2.4; 2.2 

0.6; 0.7 

0.7; 0.7 

0.1; 0.2 

0.2; 0.2 

BH-BH 
initial design 
xylophone 

29.9; 34.1 
26.6; 29.9 

69.0; 69.5 
58.8; 58.6 

25.6; 29.5 
22.9; 26.0 

2.3; 2.8 

2.1; 2.5 

TOTAL 
initial design 
xylophone 

31.7;35.9 

29.0;32.2 

76.3;77.0 

70.6;70.4 

26.8;30.6 

24.6;27.6 

3.0;3.4 

3.3;3.5 




Figure 2. The observed lensed GW event number distribution as a function of z. Left figure 
corresponds to pintr < 8 for the /_ image. The right one corresponds to the /_ image including both 
Pint < 8 and pint > 8. “Low-end” metallicity galaxy evolution and standard model of DCO formation 
are assumed. 


for the I- image to be detected. This guarantees that the /+ image will be observed, too. It 
means that in principle one would be able to establish the nature of such two time-delayed 
signals with similar temporal structure (frequency drift) but differing only in amplitudes, as 
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Table 2. Expected numbers of lensed GW events with pmtr < 8 for which the /+ image is magnihed 
above threshold po = S. Other assumptions and terminology - like iir Table 1. 


DCO scenario 
metallicity evolution 

standard 
High; Low 

optimistic CE 
High; Low 

delayed SN 
High; Low 

high BH kicks 
High; Low 

NS-NS 
initial design 
xylophone 

2.1; 1.5 

3.9; 3.2 

17.9; 20.4 
40.5; 43.5 

2.4; 1.6 

4.4; 3.5 

2.3; 1.5 

4.1; 3.2 

BH-NS 
initial design 
xylophone 

6.4;7.1 

7.5; 7.9 

11.8; 11.3 
12.8; 11.9 

3.3; 3.6 

3.7; 3.8 

0.7; 0.8 

0.9; 0.9 

BH-BH 
initial design 
xylophone 

161.8; 184.1 
144.2; 161.9 

373.4; 376.2 
318.5; 317.4 

138.3; 159.3 
124.0; 140.7 

12.5; 14.9 

11.5; 13.6 

TOTAL 
initial design 
xylophone 

170.4;192.7 

155.6;173.0 

403.0;407.9 

371.8;372.9 

144.0;164.5 

132.0;148.0 

15.5;17.2 

16.5;17.7 



Figure 3. Probability density of DCO inspiral events yearly rate (as a function of redshift) to be 
observed by the ET in its initial design. Continuous lines refer to the total catalogue of lensed and non- 
lensed systems, non-lensed systems are represented by dashed lines. It illustrates the magnification 
bias for different DCO systems. Note the logarithmic scale used in this picture. In particular figure 
shows that magnification bias is negligibly small, and even not noticeable for BH-BH systems. 


a consequence of gravitational leasing. Hence, the rates shown in Table 1 supplement the 
rates reported in Table 3 and Table 4 of Biesiada et al. [5] . Concerning the total expected 
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rates of lensed GW events, one can see that they are roughly doubled for the initial design 
of ET, for the “xylophone” design this effect is smaller (about 50 — 60% increase) but still 
substantial. Therefore the numbers displayed in Table 1 quantify magnification bias in the 
catalogue of lensed GW inspiral events. It is most pronounced for NS-NS systems, where one 
has an order of magnitude increase in the expected lensing rates. 

Figure 1 shows the distribution of yearly rates of lensed events with respect to the 
intrinsic SNR parameter p in different classes of DGO within the reach of ET detector in its 
initial and “xylophone” design. One can see that NS-NS population peaks below pintr. = 5 
with a small high-end tail of the distribution being intrinsically stronger than the detection 
threshold po = 8- This is reflected in marginal contribution of NS-NS systems to the lensing 
rates discussed in [5] and explains why these systems are affected the most by magnification 
from strong lensing. However, even for the BH-BH binary systems a noticeable part of 
the distribution lies below the threshold. Figure 2 displays the distribution of yearly rates 
of lensed events with respect to the source redshift Zg- Left panel shows the systems with 
intrinsic SNR below the threshold, while the right one shows the combined distribution of the 
total available DGO population (having SNR below as well as above the detection threshold). 
One can see that gravitationally lensed intrinsically faint sources probe higher redshifts, so 
the future catalog of gravitationally lensed GW events would be contaminated by higher 
redshift sources — in agreement with the general idea of how magnification bias works. This 
broadens the ET’s distance reach and may enable us listen to the waves of DGO inspirals 
from redshifts z = 2 - 8 which will lead to better understanding of the early epochs of star 
formation. 

One of the most important issues in modern observtional cosmology is the determination 
of star-formation history. Over last two decades a coherent picture emerged [13] according to 
which star formation rate (SFR) increase to the redshift z=2 {SFR[z = 2) ~ 10 x SFR{z = 
0)) and then decreases. This was possible to achieve because of massive spectroscopic surveys 
like SDSS and invaluable information the spectrum of a galaxy bears concerning its star 
formation history and gas metallicity. However the behavior of SFR at higher redshifts is 
still uncertain. It has been know that coalescence rates probed by GW laser interferometric 
observations could be used to test alternative scenarios of star formation. For example the 
ET Design Study document demonstrates that the ET will be able to discriminate between 
four SFR models: Hopkins & Beacom [14], Nagamine et al. [15], Fardal et al. [16] and 
Wilkins et al. [17]. The DGO coalescence detection rate would be the best source of such 
constraints, lensed events of intrinsically faint sources will not add much to it. However, 
still remains the possibility to measure the luminosity distance and redshifted chirp mass of 
such sources. Assuming some reliable background cosmology one would be able to estimate 
the intrinsic chirp mass of the system. The possibility to determine DGO masses up to z=5 
(or even higher for BH-BH systems) is intriguing. It is because such distant DGO systems 
are fossils of the era of high-mass star formation in the Universe. Let us also remind that 
detection of coalescing DGO at certain redshift yields information about BH and NS formed 
at even earlier epochs because of the delays between formation and coalescence. Measuring 
the masses of DGO containing BH will open a new chapter of astrophysics and certainly will 
tell us a lot about star formation scenarios. The comprehensive insight into the underlying 
distribution of NS masses would provide the means to study not only specific particular 
aspects like NS matter equation of state but also more fundamental ones. For example as 
the masses of NSs also retain information about the past value of the effective gravitational 
constant G, with the determination of the NS mass range at high redshifts it may be even 
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possible to probe the potential evolution of such physical constants [18]. Equally intriguing 
question is about the influence of dark matter accumulating in the cores of the NS on their 
properties (masses, conversion to quark stars etc.). This issue has usually been neglected in 
astrophysical studies, but recently attracted growing attention [19-21]. In particular Kouvaris 
& Angelez Perez-Garcia [22] try to explain the NS braking index by dark matter whereas 
Fuller & Ott [23] invoke dark matter induced collapse of NS as a link between fast radio 
bursts and missing pulsar problem. 

In addition to the magnification bias on the lensed events, one can estimate the magni¬ 
fication bias at the level of full DCO inspiral events catalogue. For this purpose one should 
calculate the detection rate of intrinsically faint events for which only I_|_ image was mag¬ 
nified above the threshold. This is shown in Table 2. From Eq.(2.7) and Eq.(2.8) one can 
see that elementary cross sections for /_|_ image magnified above threshold are higher than 
in the case of /_ image. Therefore the numbers reported in Table 2 are much higher than 
analogous numbers in Table 1. However, they should be compared with yearly detection 
rates of DCO inspiral events predicted for the FT (Table 1 and Table 2 of Biesiada et al. 
[5] ). Such comparison shows that the magnification bias at the level of the full inspiral event 
catalog would be 0.001. Therefore this would not affect much cosmological inferences drawn 
from such catalog. Figure 3 illustrates this effect by plotting together probability density 
of yearly detection rate of non-lensed (dashed line) sources and total prediction — enriched 
by systems with /_|_ image magnified (solid line). In particular, one can see how the lensing 
effect extends the high redshift tails of these distributions for all DCO systems except the 
BH-BH ones. However, the logarithmic scale was adopted in order to visualize this effect 
and the above mentioned claim of negligible magnification bias remains valid. 

4 Conclusions 

We all hope that the new era of gravitational astronomy will soon be opened by the second 
generation of (upgraded) interferometric gravitational wave detectors. The expected benefits 
of direct detection of GW from astrophysical sources — inspiralling DCO systems in the 
first place — are much more than just seeing GW in flesh. They will provide unique tests 
for fundamental physics (e.g. alternative theories of gravity) and invaluable complementary 
tests of relativistic astrophysics (stellar and supermassive black holes) and cosmology. All 
the above mentioned is even more true for the ET which is expected to provide an extremely 
rich catalog of DCO inspirals detecting up to thousands of them per year. This way, ET will 
be sensitive to a population of sources at very high redshifts, allowing to study cosmological 
evolution of sources, the history of star formation and its dependence on the matter content 
of the Universe, and the development of large-scale structure in the Universe. 

In this paper we extended our previous studies [4, 5] on gravitational lensing of such 
DCO inspiral events. In particular we explored the case when DCO systems are intrinsically 
faint, i.e. having SNR parameter p < 8. This means that usually they would not be detected 
by ET, but the presence of gravitational lenses along the line of sight changes this situation 
and they could become visible to ET being magnified by lens. Our main result is that tens 
to hundreds of such extra events could be detected per year. Especially NS-NS (but also 
NS-BH) systems are affected by this mechanism and the ET’s distance reach for signals 
from such coalescences broadens from z ~ 1 in the non-lensed case to z ~ 4 in the lensed 
case. This opens possibility to study star formation history covering substantially earlier 
epochs of the Universe. Lensed DCO events could also give us a unique information about 
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masses and NS equation of state at much higher redshifts than one would be able to test 
without magnification by lenses. This inference would be based on individual systems not 
on statistical reasoning, so any single faint DCO inspiral would be useful. However, equally 
important will be to determine the Hubble constant or cosmic equation of state using catalogs 
of inspiral events. Therefore, one may worry about the adverse effect of contaminating the 
catalogs by lensed events (so called magnification bias). Hopefully our results show that the 
magnification bias is negligible and should not affect much cosmological inferences. 
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